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Abstract 



We examine isovector and isoscalar neutron-proton correlations in an ex- 
actly solvable model based on the algebra SO (8). We look particularly closely 
at Gamow- Teller strength and double beta decay both to isolate the effects 
of the two kinds of pairing and to test two approximation schemes: the renor- 
malized neutron-proton QRPA (RQRPA) and generalized BCS theory. When 
isoscalar pairing correlations become strong enough a phase transition occurs 
and the dependence of the Gamow- Teller (5 + strength on isospin changes in 
a dramatic and unfamiliar way actually increasing as neutrons are added to 
an iV = Z core. Renormalization eliminates the well-known instabilities that 
plague the QRPA as the phase transition is approached, but only by unnat- 
urally suppressing the isoscalar correlations. Generalized BCS theory on the 
other hand, reproduces the Gamow- Teller strength more accurately in the 
isoscalar phase than in the usual isovector phase, even though its predictions 
for energies are equally good everywhere. It also mixes T = and T = 1 
pairing, but only on the isoscalar side of the phase transition. 

21.60.Fw, 23.40.-s, 23.40.Hc 
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I. INTRODUCTION 



Pairing correlations are an important feature of nuclear structure [rj. In heavy nuclei 
such correlations between neutrons and protons are usually neglected on the grounds that 
the two fermi levels are far apart. In nuclei with N ~ Z, however, the fermi levels are close 
and neutron-proton (np) pairing correlations can be expected to play a significant role in 
nuclear structure and decay (for a review of work on np pairing theory see Ref. f2|). The 
importance of these proton-rich nuclei in astrophysical nucleosynthesis makes it vital that 
the np correlations are well understood, and upcoming experiments with radioactive beams 
will soon test our understanding. 

Taking np correlations seriously complicates the usual treatment of pairing, which 
stresses the interaction of like particles in time-reversed orbits, i.e. the formation of pp 
and nn pairs. Generalizing this picture raises at least two new issues. First, the pp, nn, 
and np isovector pairs must all be treated on an equal footing so that isospin symmetry is 
respected as much as possible 0. Second, the competition between two kinds of np pair- 
ing — isovector and isoscalar (T = 1 and T = 0) — must be taken into account. This 
issue apparently arises even in nuclei with N > Z where np pairing is by most measures 
small. For example, the rate of two-neutrino double-beta decay within the np Quasiparticle 
Random Phase Approximation (QRPA) is extremely sensitive to the strength of isoscalar 
particle-particle (i.e. pairing) correlations, making reliable calculations difficult. When these 
correlations become strong enough the method fails even to give finite answers. 

The breakdown in the QRPA signals an impending phase transition. Is it real or an 
artifact of the assumption underlying the approximation that the ground state contains no 
np correlations? What are the properties of the "isoscalar phase" , if it is real? Is the Renor- 
malized QRPA (RQRPA) |4]||, in which solutions are more stable, a good way to handle the 
breakdown/phase transition? Is generalized BCS theory, a scheme for treating np pairing 
on a more equal footing with nn and pp pairing, able to quantify the interplay between 
the two phases? We address these questions here in a solvable model that incorporates 
both isovector and isoscalar pairing, making it considerably richer than the more schematic 
models (e.g. the Lipkin model JJ) typically used for this kind of study. 

The structure of this paper is as follows: In section 2 we describe the model and its 
analytic solution for energies and Gamow- Teller beta-decay matrix elements, stressing the 
existence of two limiting solutions corresponding to pure isovector or isoscalar pairing. We 
show that in the isoscalar phase, charge-changing processes have counterintuitive features. 
Section 3 contains an outline of the QRPA and RQRPA as realized in the solvable model, 
and applies them to single and double beta-decay to test the quality of the approximations. 
In Section 4 we describe generalized BCS theory for 50(8), again with emphasisis on beta 
decay strengths (though we also examine ground-state energies), and again test the reliability 
of the approximation scheme. Section 5 is a conclusion. 

II. THE MODEL AND ITS EXACT SOLUTION 

We consider a set of degenerate single particle orbitals, characterized by I, s = 1/2, t = 
1/2. The total number of single-particle states is Q = J2i{2l + 1)- We make the model 
solvable by building a basis entirely from L = operators: 5^, which creates pairs with spin 
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5 = and isospin T = 1 (with projection is), and Pj, which creates a pairs with 5* = 1 and 
T = (spin projection //). Together with the one body operators that generate SU(4) - 
the total spin 5, the isospin T, and the operator JF^ = a{i)^r{i) v - the 7 = pair 
creation and annihilation operators form the algebra 50(8) J7|- The physics associated with 
this model has been studied previously [0-0, but with emphasis on energy levels; here our 
focus will include charge-changing decay 

The most general Hamiltonian invariant under 50(8), omitting terms such as 5 • 5 and 
T • T that affect energies but not wave functions, depends on three parameters and has the 
form H 

h = E sts v - 9 -^± E pIp, + g#rP H • (i) 

The first term in the Hamiltonian corresponds to isovector spin-0 pairing, the second rep- 
resents isoscalar spin-1 pairing, and the last is a (primarily) particle- hole force in the T = 1 
5=1 channel. 

In certain important limits, analytic expressions for energies and wave functions have 
been derived. If x — 1 and g p h = 0, the Hamiltonian @ conserves an 50(5) subalgebra 
and corresponds to "standard" spin-singlet isovector pairing, with np pairs treated on an 
equal footing with like-particle pairs [|K],[iT|]. The eigenstates, characterized by the number 



of nucleon pairs M (we consider only nuclei with an even number of nucleons), the isospin 
T, and the singlet-pairing seniority v s , have energies 

E{v a , T) = -f [(2JV - v s )(4{l + 6 - 27V - v s ) - 4T(T + 1)] . (2) 

o 

Similarly, for x — — 1 and g p h = the exact solutions are characterized by the spin 5 and 
the triplet-pairing seniority v t , and an analogous formula applies with T — > S,v s — > v t . 
This is the "isoscalar phase" that will cause the breakdown of the QRPA. Finally, if x = 
the Hamiltonian is invariant under 577(4). The eigenstates are then labeled by a quantum 
number A corresponding to the irreducible 577(4) representation [A, A, 0] as well as by 5 and 
T, and the eigenvalues are 

E{\, 5, T) = -| [2Af(n + 3) - M 2 - A(A + 4)] + g ph [A(A + 4) - 5(5 + 1) - T(T + 1)] . 

(3) 

5 + T must be even if TV is even and odd otherwise (this is true no matter what the 
Hamiltonian). The quantum number A has values A = 5 + T, 5 + T + 2,..., \mox> where 
^Max — TV if jV < and 2Q — Af otherwise. 

The eigenvalues and eigenstates of the general hamiltonian in Eq. (1) can be obtained 
by diagonalizing in the 5i7(4) basis. The matrices are tridiagonal and have very small 
dimension. Expressions for the matrix elements (with several typos) appear in Ref. |§, 
which, however, ignores the particle-hole interaction. The same model with the particle- 



hole interaction included was solved approximately in Ref. [ |12jj . 

By varying the parameter x one can study the phase transition from the standard spin- 
singlet isovector pairing phase, through the SU(4) Wigner supermultiplet phase, into the 
spin-triplet isoscalar pairing phase. In Fig. 1 we show the effects of this transition (with g p h 
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fixed at zero) on the overlap between the ground state and the ground state of the standard 
spin-singlet paired system with x = 1 and g p h = 0. [The abcissa is labeled by —x, so 
that the more familiar isovector phase is on the left.] The change produced by finite g p h is 
illustrated in Fig. 2; when g p h > increases, the Hamiltonian more nearly conserves £77(4) 
symmetry and the phase transition becomes less pronounced. The overlap is an obvious 
"order parameter" in the model, and its point of inflection locates the phase transition. In 
even-even systems (£ = 0, T even) this point shifts to the right from the ££7(4)-limit value 
x = as T increases. The reason is that any excess neutrons are necessarily in isovector 
pairs, making the transition to isoscalar pairing more difficult for the remaining nucleons. 

We turn our attention now to transitions induced by the GT operators Their matrix 
elements are easily evaluated when the wave function is written in the basis \J\f, A, S, T) 
because the operator is diagonal in Af and A and can change £ and T by one unit only. 
Explicit formulae can be constructed from the SU(4)/SO(4) Clebsch-Gordan coefficients 
derived in Ref. flT3"f . The GT operators either increase or decrease N — Z (assumed to be 
nonnegative); the corresponding strengths are called [3 + and (3~ . The two strengths are 
constrained by the Ikeda sum rule 

£(/T ) - S((3 + ) = 3(N - Z) . (4) 

In the S77(4) limit the (3 + strength vanishes and one state exhausts the f3~ strength. In the 
two extreme limits surrounding £77(4), i.e. x — 1, g p h = (the isovector pairing phase) and 
x = — 1, g p h = (the isoscalar phase) we can derive analytic expressions for £(/3 + ) from the 
ground state as a function of T = T z = (N — Z)/2, since in those limits the Hamiltonian 
contains only the generators of an £0(5) subgroup. The /5-decay operators break one pair 
(of either kind), leading to matrix elements between simple £0(5) representations, the 
properties of which were studied in Ref. fl4]| . In the spin-singlet isovector phase we find 

c (f}+) = W-T)(T + l)(2n-AT-T) 
^ ' (2T/3 + l)(fi + l/2) ' y ' 

This result applies in a single j-shell with degeneracy 2Vt as well as in the model discussed 
here, which necessarily contains at least 2 degenerate levels in the j—j scheme (j = /±l/2). 
The (3 + strength is plotted as a function of T for Q = 12, M = 10 in Fig. 3. Except for the 
initial plateau at low T the behavior of the curve is qualitatively similar to that obtained 
in BCS theory, where the neutron-proton interaction is ignored. The gradual decrease in 
strength with T is caused by Pauli blocking. 
In the isoscalar phase the result is 

S{P } = Q-T+l/2 • (6) 

Here the behavior of the strength as T increases from is surprisingly different (see Fig. 3). 
The substantial rise at first seems counterintuitive since the neutron excess is increasing. 
Blocking is not the only factor at work, however. For lower T the effect is overcome by the 
collective behavior of the boson-like £ pairs in the final state. 

In Fig. 4 we examine the behavior of the strength S((3 + ) between the two limits, as a 
function of — x for fixed g p h = (When g p h ^ all curves become flatter, because the system 
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is closer to the 577(4) limit.) As — x increases, the strength S(f3 + ) decreases, vanishing when 
the S77(4) limit is reached and increasing again as the isoscalar pairing phase is approached 
until finally it is considerably larger than in the isovector phase. The large strength is 
caused in part by the transfer of protons from pp pairs, which cannot participate in (3 + 
decay, to isoscalar np pairs, which can. Only close to x = —1, however are this effect and 
the parabolic isospin dependence fully present; when — x is small the strength can be small 
as well and the isospin dependence complicated, as is apparent from the crossings of curves 
in Fig. 4. A large (3 + strength (compared, e.g., to the Ikeda sum rule) therefore reflects 
very strong isoscalar pairing. If real, it would have important consequences for r-p process 
nucleosynthesis. 



III. DOUBLE-BETA DECAY, THE QRPA, AND THE RQRPA 

We have stressed (3 + strength because of its simplicity and sensitivity to details of nuclear 
structure. Now, however, we want to discuss modifications of the QRPA, the most frequent 
and controversial application of which is to double-beta decay. Actually, in a model with as 
few states as this one, the (3 + strength from the "final" nucleus (/) essentially determines 
the double-beta decay matrix element (considered here in the closure approximation for 
simplicity and because the energy denominator can change without a concommitant change 
in the wave functions). The reason is that the matrix element has the form 

M 2 J T (d) = Wl £<?(*) •<?0>(z)_T(j)_|0+> • (7) 

and for a moderate neutron excess the f3~ strength, the other relevant quantity, hardly 
varies with x. In realistic calculations, the QRPA, which has many desirable features, 
suffers an unfortunate instability when g pp = (1 — x)g/2 becomes too large that manifests 
itself through infinite values for both the (3 + strength and the double-beta-decay rate. A 
number of remedies have been proposed recently. One that has received particular attention 
is the Renormalized QRPA (RQRPA) which eliminates the instability of the QRPA 

through a self-consistent calculation of the ground state. The model presented here is ideal 
for examining how both QRPA and the new approximation work. 

The p-n QRPA, described, for example, in Ref. [0 and applied to SO (8) in the same 



paper, begins with the ordinary BCS ansatz, a coherent state of isovector neutron-neutron 
and proton-proton pairs, and proceeds by admixing neutron-proton quasiparticle pairs into 
the ground and excited states. In SO (8) the procedure leads to 2 by 2 matrix equations in 
each of the two (Fermi/isovector pairing and Gamow-Teller/isoscalar pairing) channels and 
can be solved by simple diagonalization. More specificially, the one excited state in each 
channel (S = or S = 1) is written in the form 

\S >= (X s [alal] L =°> s - Y s [a p a n } L=0 ' s ) |0 > , (8) 

where a^ p , ot^ (a p , a n ) create (destroy) proton and neutron quasiparticles, the brackets 
indicate angular momentum coupling, |0 > is the QRPA ground state and (X ,Y ), and 
(Xi,Yi) are the "physical" eigenvectors in the spin-0, isospin-1 (Fermi) and spin-1 isospin-0 



(Gamow- Teller) channels (additional details are in Ref. ||12|| ). The two channels decouple 
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and for two-neutrino double beta decay only the second is relevant. Associated with the 
S = 1 eigenvector is an eigenvalue that becomes complex when g pp reaches a critical value 
connected with the impending phase transition. The states with complex eigenvalues are 
not normalizable and have no physical significance, so that the approximation fails to give 
even an incorrect answer beyond the critical point. This is the "collapse" referred to above 
and is preceded by rapid changes in the f3 + and double-beta-decay amplitudes. 

In the RQRPA, described for charge-changing modes in Ref. the two channels are 
coupled in an attempt to make the vacuum self consistent, and the resulting equations are 
nonlinear. For the 5*0(8) model the equations have seven variables: the two sets of X's and 
F's, the eigenvalue associated with each set, and a renormalization parameter. The iterative 
procedure advocated in Ref. {§ often does not converge here, but the model's simplicity 
makes the equations easy to solve by other means. Unlike the QRPA, the RQRPA never 
exhibits the analog of the complex eigenvalues that signal instability, and therefore never 
yields rapidly changing matrix elements. The question is whether any important physics is 
lost in the process of guaranteeing a ground state that is built on the BCS state. 

Fig. 5 presents the exact (3~ and f3 + strengths for fixed M and Q and several values 
of T, along with the QRPA and RQRPA approximations to the strengths. The QRPA 
breakdown is reflected in the expression for the f3 + strength, which blows up at the critical 
value of x (or g pp ). By contrast the RQRPA strength is perfectly stable. The graphs 
make it clear, however, that the stability is achieved at a significant price; the very real 
phase transition to a ground state dominated by isoscalar pairing correlations changes the 
behavior of S(j3 + ), causing the QRPA to break down, but refuses to show itself at all in the 
RQRPA approximation. The reason is that in preserving (self-consistently) the basic QRPA 
ansatz the RQRPA limits the isoscalar correlations in the ground state. Put another way, 
the QRPA breaks down for a reason; there really is a phase transition and it really is nearby, 
and the RQRPA erases all traces of it. Thus at the very point at which the QRPA fails the 
RQRPA also begins to deviate badly from the exact result. To make matters worse, and 
this has been noted elsewhere |H| , renormalization destroys one of the nicest features of the 



QRPA, the preservation of the Ikeda sum rule Eq. (4). In this model, at least, nothing is 
gained by using the RQRPA. 

To demonstrate this explicitly for double beta decay, we show in Fig. 6 the matrix element 
Mq u t for Af = 12 and T = 4 as a function of g pp /g P air, where g pa i r = (1 + x)/2. We use this 
parameter rather than x because it more closely resembles that used in realistic calculations. 
We have set g p h = 1.5g so that the QRPA breaks down just beyond the point at which the 
matrix element crosses the origin (when g pp = g pair i.e. at the S77(4) point). This is the 
situation in more realistic calculations as well, but our model shows it to be pure coincidence; 
the breakdown of the QRPA moves to larger g pp as the essentially independent parameter 
g p h is increased, while the crossing point never moves, implying that nothing fundamental 
is behind the proximity of the crossing to the point at which the QRPA fails in realistic 
calculations. Interestingly, the exact matrix element in our model varies smoothly as the 
phase transition is traversed, in sharp contrast to the rapid drop in the prediction of the 
QRPA, which as noted blows up completely just past the crossing point. The RQRPA, as 
in the last figure, begins to fail at the same point and offers little advantage over the QRPA 
itself. 

In more realistic calculations this last conclusion may or may not hold. The simple model 
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examined here contains only two very collective degrees of freedom. It is certainly possible 
that with less collectivity the QRPA approximation is worse and the breakdown occurs 
far from the actual phase transition. In that event the RQRPA would offer advantages, 
especially in the region between the breakdown of the QRPA and the real phase transition. 
It would therefore be useful to examine the approximation in a model that dilutes the 
collectivity of the T = np pairs but is still solvable. A two- level version of SO(8), i.e. 
50(8) x 50(8) P|, might be a good place to start. Here we can say only that we find no 
evidence supporting the validity of the RQRPA. 



IV. APPLICATION OF "GENERALIZED BCS THEORY" 

In this section, we apply generalized pairing theory to the 50(8) model, to assess its 
ability to provide a meaningful approximate description of the ground-state dynamics of 
the model in the various phases. We simplify the model slightly by setting g p h to zero (and 
g to 2, which merely scales the energies). 

Generalized pairing theory is well reviewed in Ref. @ and thus will not be discussed in 
detail here. Suffice it to say that the theory is founded in the HFB approximation, sup- 
plemented by the further assumption that the only nonzero matrix elements of the Hartree 
Fock and pair potentials are those connecting the four states \ap >, \an >, \ap > and 
\an >, where \a > denotes the state obtained by time-reversal on the state \a >. As such, 
the theory naturally accomodates all pairing modes on an equal footing. This includes the 
usual pp and nn pairing as well as pn, np and pn pairing. Note that here and in subsequent 
discussion we explicitly distinguish the pairing of particles in "the same orbit" (e.g., np) 
from the pairing of particles in "time-reversed orbits" (e.g., np). 

We have chosen to formulate the theory in terms of the density matrix p(a) and 
the pairing tensor t(a). These matrices, after invoking time-reversal invariance, take the 
parametrized forms 
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with the coefficients interrelated by four unitarity conditions 

(1 - pi - p 2 )po - (h + t 2 )t = , 



Pi 



i - pi - p 2 )p 3 + (h + h)t 3 = o , 
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(10) 
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In our application to the 50(8) model, we impose constraints on the average number of 
neutrons and the average number of protons of the system, thereby fixing the parameters 
px and p 2 according to 
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Two more parameters, po and p$, are fixed from the first two of the unitarity conditions 
(0)- 

Our procedure is first to express the energy of the generalized quasiparticle vacuum as a 
function of the remaining five parameters of the density matrix and pairing tensor and then 
to look for local minima, rather than to solve the usual self-consistent eigenvalue equation. 
The remaining two unitarity conditions (|10D are implemented via Lagrange multipliers. The 
system of equations arising from these variational conditions in principle admits several 
solutions (all having 9 = or vr/2), the energetically- lowest of which defines the generalized 
BCS approximation to the ground state of the system. 

Fig. 7 shows the energies associated with the solutions to the generalized BCS equations 
for the case of Q = 12 and M = 5. The results are plotted as a function of the hamiltonian 
parameter x and for various values of the neutron number N. The solutions displayed in 
the figure have the following character: 

• Solution A corresponds to pure pp and nn pairing. 

• Solution B corresponds to pure T = 1 pn and np pairing. It only exists when N=Z, 
where it is precisely degenerate (but does not mix) with solution A. 

• Solution C in general involves both T=0 pn and np pairing and T=l pp and nn pairing. 
The relative importance of these different pair correlations is dictated by the three 
parameters to, t\ and t%. To a good approximation, the first reflects the number of 
collective T=0 pn and np pairs, whereas the latter two reflect the number of collective 
pp and nn pairs, respectively. As x increases from -1, the solution eventually merges 
into solution A, ceasing to exist beyond that "critical point". At precisely this point, 
there is a change in the character of the ground state predicted in generalized BCS 
approximation that mirrors the true state of affairs. 

Figure 8 shows the generalized pairing results for the ground state energy in comparison 
with the exact energies discussed earlier, again for Q = 12 and Af = 5. Included are 
results corresponding to values of the neutron particle number N ranging from N = 5 to 
N = 10. The results for N = through N = A follow from the symmetry of the problem. 
The generalized pairing results correctly reproduce the trends of the exact results, equally 
well in both phases. The good predictions hold up for even-even nuclei, odd-odd nuclei, 
symmetric nuclei with N = Z, and nuclei with N = or Z = 0. The generalized pairing 
approximation even reproduces the gradual shift of the phase transition to negative values 
of x when the difference between N and Z increases. 

We now turn one last time to Gamow- Teller matrix elements, focusing here (as above) on 
the summed f3 + strength. We obtain the strength by evaluating the quasi-particle vacuum 
expectation value of the 1+2 body operator 

S/3+ =3 J2l 1 , mil , msi ,l 2 , mi2 , ms2 ,ii {~ Y (1 A* 2 m sil| m si — A 4 2 m s2\\ m s2 ~ A 4 ) 

Pl2,mi 2 ,m S2 —iJ, n h,mi 2 ,m 32 n l lt m llt m si +fi Ph,mi 1 ,m ai • (12) 
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Here p^ (p) creates (annihilates) a real proton and n) (n) creates (annihilates) a real neutron. 
When the ground state is dominated by ordinary pp and nn pairing and represented by 
solution A, the total (3 + strength is given by the usual formula 

V = SZ - ?™ . (!3) 
When solution C applies, the result is 

V = 3Z - — + 4fi 2 po - 2fi(* ) 2 • ( 14 ) 

(The (3~ strengths follow from the above and the Ikeda sum rule, which is preserved at the 
BCS level). In Figure 9, we compare the exact (Eqs. (|5|) and (g)) and generalized BCS results 
for the summed (3 + strength as a function of T z in the isovector and isoscalar limits when 
Q = 12 and M = 10. Solution A, which is equivalent to the ordinary BCS product wave 
function, reproduces the general trends of the exact (3 + strength but differs significantly as 
T z — > 0. The results in the isoscalar limit are much better for all values of T z , including 
T z = 0. Only for energies are the BCS results equally good in both phases. 

Figure 10 compares the exact and generalized BCS f3 + strengths for the same Q and 
M values as in Figure 9. Here, however, we only consider the results for N = 14 and 
Z = 6 [T z = 4), but for all values of x. Once again, the generalized BCS approximation 
provides a somewhat better reproduction of the trends in the exact results when the system 
is dominated by isoscalar pairing correlations. The suppression of f3 + strength that shows 
up as the SU (4) limit is approached, however, is not described well by the generalized BCS 
formalism, either on the isoscalar or isovector side. More accurate results in this regime may 
require an extension beyond our BCS approximation that fully accomodates the coexistence 
of all the different pairing modes on both sides of the phase transition. 

A striking feature of Fig. 10 is the rise of f3 + strength in the isoscalar region as — x 
increases to 1. That this rise is correlated as noted earlier with the transfer of nucleons from 
pp and nn pairs to T = np pairs (in time reversed orbits — we have dropped the overscore) 
is apparent from Fig. 11 where we plot the number of collective pairs of different types for 
the same system as in Fig. 10. We show both the generalized BCS results and the "exact" 
results, using the standard prescription for operators that roughly measure the number of 
collective pairs (It is this prescription, which we believe can be improved by a better 
treatment of Pauli effects, that is responsible for the appearance of low-levels of np pairs 
before the phase transition, where the corresponding pairing tensor is zero.) Part b) of the 
figure clearly shows that in BCS theory the wave function does not change from its product 
form until the phase transition is reached. It also shows that the constant wave function 
is actually not so bad an approximation when T z is as large as in the figure. Only when 
iV m Z does the theory, which will not reflect np pairing to the left of the phase transition, 
fail badly. Even then, however, the BCS approximation manages to reproduce energies well. 



V. SUMMARY AND CONCLUDING REMARKS 

We have examined the interplay between isovector and isoscalar pairing modes in an 
exactly solvable SO(8) model, focusing on the matrix elements of charge-changing opera- 
tors. The behavior of the fl + strength in the isoscalar phase is counterintuitive, rising with 
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increasing neutron excess instead of falling. Double-beta decay, by contrast, varies smoothly 
and predictably on both sides of the phase transition. This behavior is in sharp contrast to 
the predictions of the QRPA, which fails when isoscalar pairing becomes too strong. 

Partly for this reason, we have tested two approximation schemes that purport to better 
accomodate neutron-proton correlations. One, the RQRPA, works well in the isovector phase 
but fails completely to capture the physics of the phase transition. Ironically the second 
approximation, provided by generalized BCS theory, does a better job for the total beta- 
decay strength in the isoscalar phase than in the isovector phase. (It successfully reproduces 
ground state energies everywhere, however, even in the vicinity of the transition). The reason 
is that the strength operator is a scalar in space and spin, and is therefore not sensitive to 
the spin "deformation" that inhabits the BCS wave function in the isoscalar phase. On the 
other hand, the isospin deformation in the isovector phase distorts the expectation value 
of the strength operator, which contains isoscalar, isovector, and isotensor pieces. In Ref. 
H it was shown in a simpler model, based on 5*0(5) and containing only isovector pairing 
interactions, that projection of the generalized BCS quasiparticle vacuum onto states with 
good isospin after variation can fix this problem (though there the analog of solution B was 
used as the "isointrinsic state"). It is far from clear, however, that projection will allow the 
dynamical mixing of isoscalar and isovector pairing before the phase transition is reached. 
Recently |[n|, the Lipkin-Nogami method was shown to do the trick at least in part, and it 
would be interesting to test it in a model like this one. On the other hand, the phases are 
mixed even at the BCS level on the isoscalar side of the critical point. 

The shortcomings of the RQRPA and the successes of generalized pairing theory raise 
the following question: Can the generalized quasiparticle vacuum be used as a starting point 
for a "generalized QRPA" that works even in the region of the phase transition? Something 
along these lines has been attempted in Ref. |T7[, but only after forcing the BCS to mix 



isoscalar and isovector pairing in an artifical way. Our BCS also appears not to mix the two 
kinds of pairs except to the right of the critical point, a region that is probably unphysical 
in nuclei that undergo double-beta decay. Perhaps a more fruitful approach, therefore, will 
be a more self-consistent QRPA, in which the RPA and BCS equations are coupled. Such 



a procedure can rescue the Ikeda sum rule [18| and could conceivably facilitate isoscalar- 



isovector mixing at the BCS level even to the left of the critical point. Other modifications of 
the basic BCS procedure, including approximations to projection, may also mix the pairing 
modes without complicating the method too severely. Which of these approaches is the 
simplest and most useful remains to be seen. 
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FIGURES 



FIG. 1. Overlaps between the ground state and the pure isovector spin-singlet paired state vs. 
the parameter x in the hamiltonian of Eq. (1), for Q = 12,7V = 10 and S = T = (solid line), 
S = 0, T = 2 (long dashes), and S = 0, T = 4 (short dashes). Here and in Figs. 2, 4, 5, and 10 the 
quantity — x is used on the abscissa axis so that the standard isovector phase is on the left. 

FIG. 2. Same as Fig. 1, with S = T = 0, except for different values of g p h- The solid line 
corresponds to g p h = 0, the long dashes to g p h = l-0g, and the short dashes to g p h = 2.0g. 

FIG. 3. The Gamow-Teller /?+ strength B(GT) vs. the initial neutron excess T z = T for 
VL = 12,7V = 10. The solid curve is for the pure isoscalar spin-triplet pairing phase and the dashed 
curve for the standard isosvector spin-singlet phase. 

FIG. 4. The Gamow-Teller (3 + strength B(GT) vs. the hamiltonian parameter — x, for 
SI = 12, TV = 10, S = and g p h = 0. The solid curve corresponds to T = 0, the long dashes 
to T = 2, and the short dashes to T = 4. 

FIG. 5. The Gamow-Teller strength B(GT) for (3~ (upper parts of each panel) and /3 + (lower 
parts of each panel) vs. —x. The exact results are denoted by the solid lines and the RQRPA 
results by the dotted lines. Calculated for Vt = 12,7V = 10, S = 0, and g p h = g and several values 
of the isospin T labelling the corresponding panels. 

FIG. 6. Double beta decay matrix element M"q T vs. g pp /g pa ir- The exact results are denoted 
by the solid line, the RQRPA results by the dashed line, and the QRPA results by the dotted line. 
Calculated for SI = 12, TV = 12, S = 0, T = 4 and g ph = l.hg. 

FIG. 7. Energies associated with the different variational solutions to the generalized BCS 
equations, for O = 12, TV = 5. The panels are labelled by the neutron number N. Here and in Fig. 
8 the abscissa is the hamiltonian parameter x running from 1 to -1. The standard isovector paired 
state is on the left as in the other figures. 

FIG. 8. Comparison of exact ground state energy (solid curve) with the approximate result 
obtained in the generalized BCS approximation (dashed curve). The results are shown for = 12 
and TV = 5, and for different values of the neutron number N. 
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FIG. 9. The Gamow-Teller (3 + strength B(GT) for VL = 12 and M = 10 as a function of the 
isospin T in the two limiting phases. Solid lines represent the exact solutions and dashed lines the 
generalized BCS solutions in the isovector (bottom two lines) and isoscalar (top two lines) phases. 

FIG. 10. Exact (solid) and generalized BCS (dashed) Gamow-Teller /?+ strength B(GT) vs. 
the hamiltonian parameter — x for £1 = 12, J\f = 10 and T z = 4. 

FIG. 11. "Numbers" of different types of pairs in the exact (a) and BCS (b) solutions as a 
function of — x, for CI = 12, J\f = 10 and T z = 4. The solid line measures the number of nn pairs, 
the dashed line the number of pp pairs, the dot-dashed line the number of T = np pairs, and the 
dotted line the number (nonzero only because of the rough definition of "number" ) of T = 1 np 
pairs. 
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